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! We consider the possibility that a right-handed sneutrino can serve as the 
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source of energy density perturbations leading to structure formation in cos- 
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. mology. The cosmological evolution of a coherently oscillating condensate of 

' right-handed sneutrinos is studied for the case where reheating after inflation is 

due to perturbative inflaton decays. For the case of Dirac neutrinos, it is shown 



that some suppression of Planck scale-suppressed corrections to the right-handed 
neutrino superpotential is necessary in order to have sufficiently late decay of the 
^ I right-handed sneutrinos. cH^ corrections to the sneutrino mass squared term 

^ must also be suppressed during inflation (|c| ^ 0.1), in which case, depending 

on the magnitude of |c| during inflation, a significantly blue (if c > 0) or red 
(if c < 0) perturbation spectrum is possible. R-parity must also be broken in 
order to ensure that the Universe is not overdosed by LSPs from the late decay 
(at temperatures 1 — 10 MeV) of the right-handed sneutrino condensate. The 
resulting expansion rate during inflation can be significantly smaller than in con- 
ventional super symmetric inflation models (as low as 10^ GeV is possible). For 
the case of Majorana neutrinos, a more severe suppression of Planck-suppressed 
superpotential corrections is required. In addition, the Majorana sneutrino con- 
densate is likely to be thermalised before it can dominate the energy density, 
which would exclude the Majorana right-handed sneutrino as a curvaton. 
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1 Introduction 



The observation of neutrino masses and mass splittings, via solar and atmospheric 
neutrinos [1], strongly suggests the existence of right-handed (r.h.) neutrinos. In 
extensions of the Minimal Supersymmetric Standard Model (MSSM) [2] which can 
accomodate neutrino masses we therefore expect to have right-handed sneutrinos. In 
the cosmology of the MSSM and its extensions, Bose condensates of scalar fields such 
as squarks and sleptons form naturally after inflation [4]. These may have important 
consequences for cosmology; for example, they allow for the possibility of baryogenesis 
and leptogenesis via the Affleck-Dine mechanism [4, 5]. Since right-handed sneutrinos 
may also form condensates, it is important to consider in some detail the cosmological 
evolution of a r.h. sneutrino condensate. 

Due to the weak couphng of the r.h. sneutrinos to the MSSM fields, a condensate 
of r.h. sneutrinos will be long-lived and so may come to dominate the energy density of 
the Universe before it decays. The question of whether the r.h. sneutrino condensate 
can dominate the energy density of the Universe when it decays has recently acquired 
some importance. It has been noted [6, 7] that if a scalar dominates the energy 
density when it decays, and if that scalar is effectively massless during inflation, then 
quantum fluctuations of the scalar during inflation can in principle account for the 
primordial energy density perturbations leading to structure formation [6, 7, 8]. This 
has been labelled the curvaton scenario [9, 10]. Should the curvaton be able to account 
for the density perturbations, the parameters of the inflation model would be more 
weakly constrained than in the conventional case where density perturbations arise 
from quantum fluctuations of the inflaton. 

Thus if there exists a natural curvaton candidate such as the r.h. sneutrino, it is 
important to confirm or exclude that candidate as a curvaton^. The main goal of this 
paper is to investigate the possibility of a r.h. sneutrino curvaton. 

The masses of the r.h. sneutrinos and their coupling to the MSSM fields are 

determined by the model of neutrino masses, in particular whether they have Dirac 
""^ Models have recently been proposed where the curvaton corresponds to an MSSM flat direction 
scalar [11] and an MSSM Higgs scalar [12]. 
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masses or Majorana masses via a see-saw mechanism [1]. The Yukawa couphng of 
the r.h. sneutrino to the MSSM fields plays a fundamental role in determining the 
evolution of the r.h. sneutrino condensate, in particular its decay temperature, its 
effective mass from interacting with the background of inflaton decay products [13] 
and its rate of thermalisation/scattering from the background. 

Throughout this paper we will consider the simplest model for infiation and reheat- 
ing, coresponding to a constant expansion rate during infiation followed by formation 
of a coherently oscillating infiaton condensate and reheating due to perturbative in- 
fiaton decays. 

The paper is organised as follows. In Section 2 we discuss the cosmological enviro- 
ment due to perturbative infiaton decays. In Section 3 we consider the evolution of the 
r.h. sneutrino condensate in this enviroment. In Section 4 we present our conclusions. 

2 Cosmological Background from Perturbative In- 
flaton Decays 

After infiation we consider the infiaton S to have a mass ms and to be coherently 
oscillating about the minimum of its potential. The infiaton is assumed to decay into 
pairs of relativistic MSSM particles with initial energy ms/2. If the initial energy of 
the decay products is sufficiently large, the initial scattering rate of the inflaton decay 
products may be small enough that thermalisation only occurs once the Universe has 
expanded sufficiently for the decay product scattering rate F^c to exceed the expansion 
rate H. There are therefore two possibihties: 

(a) Instantaneous Thermalisation. The infiaton decay products thermalise imme- 
diately after decay. In this case there will be two epochs: (i) Inflaton Matter Domina- 
tion (IMD), where the energy density of the Universe is dominated by the coherently 
oscillating infiaton, and (ii) Radiation Domination (RD), defined to mean domina- 
tion of the energy density of the Universe by relativistic infiaton decay products, not 
necessarily thermalised. 
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(b) Non-instantaneous Thermalisation. The relativistic inflaton decay prod- 
ucts are unthermalised initially. In this case we will show that thermalisation cannot 
occur during IMD and so must occur during RD. Therefore in this case the back- 
ground will have three distinct epochs: (i) IMD, (ii) RD pre-thermalisation and (iii) 
RD post-thermalisation. 

2.1 Inflaton decay product thermalisation 

Prior to thermalisation there is a spectrum of decay products, ranging from red-shifted 
products from the earliest decays (occuring at the end of inflation) to products from 
the most recent decays. Most of the energy density in decay products will come from 
the most recent decays. Following the discussion of [14], the spectrum of unthermalised 
decay products as a function of energy during IMD and RD epochs is given by 

3 f Hy ps {Hr) E^/^ 
dE ~ 2 \hJ ms mT ' ^ ^ 

where 7=1 during IMD (3/4 during RD), n{E) is the number density of decay 
products with energy less than E and Hr is the expansion rate at the onset of radiation 
domination. Thus n{E) oc E^/'^ and so 80% of the decay products at a given time have 
energy between E^ and -Ed/3, where E^ is the energy of the most recent decay products. 
In addition, once thermalisation by scattering begins, the lower energy decay products 
tend flrst to increase their energy by scattering from the more numerous higher-energy 
decay products in the spectrum [15, 16], so that they may be regarded as higher- 
energy decay products as far as thermalisation is concerned. Thus we will consider 
the energy of the decay products at a given time to be approximately Ed- During 
IMD, Ed ~ ms/'2. Once the inflaton condensate has decayed away and the Universe 
enters the RD epoch, the energy of the dominant unthermalised decay products will 
be red-shifted to E^ where aR is the scale factor at the onset of radiation 

domination. 

The centre of mass (CM) cross-section of the relativistic inflaton decay products is 



where Ecm ~ 2£'d and age — 5'^/47r, where is a typical MSSM gauge or Yukawa 
couphng. Therefore the scattering rate of the infiaton decay products is Tsc = nasc, 
where n is the number density of infiaton decay products. 

During IMD most of the decay products at a given scale factor are produced during 
the previous e-folding. The energy density in the infiaton condensate during IMD is 

3 



Ps-[^) e-'^'Pi, (3) 

where r,^ is the decay rate of the infiaton, Oe is the scale factor at the end of inflation 
and Pi is the energy density during inflation, assumed constant. Therefore the number 
of inflatons which decay during an e-folding {6t ~ H~^) (and so the number of infiaton 
decay products produced) is 

(4) 

H ms 



for Tfi/H < 1. The condition for thermalisation to occur during IMD, Fgc ^ H, is then 

1, (5) 



where M — Mpi/s/Stt and we have used — ms/2 and ps ~ SM'^H'^ during IMD. 
Eq. (5) is independent of the scale factor, so if it is not satisfied immediately at the end 
of inflation, it will not be satisfied during IMD. Therefore for the case of perturbative 
infiaton decay, thermalisation must either be instantaneous or must occur during RD. 
After IMD, the number density and energy of the relativistic decay products is 



The infiaton decay rate, T^, may be expressed in terms of the reheating temperature^, 
Tr, as Fa = hrJ^^Mpi, where kr^ = {47r^g{TR)/45y/^ and ^(T^) is the effective 
number of massless degrees of freedom in thermal equilibrium [20]. (In the following 
^The reheating temperature is defined in the following to be the temperature of a thermalised 
Universe at the onset of radiation domination. It is therefore used generally to parameterise the 
energy density at the onset of radiation domination, even if the relativistic decay products at that 
time have not yet thermalised. 
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we will consider ~ 20, corresponding to the field content of the MSSM with 
g{Tn) ^ 200.) 

The thermalisation condition na,,e ^ H (where H = (^^) H{aR) during RD) then 
implies that thermalisation occurs at scale factor ath given by 



ath /":s' 



(7) 



with the corresponding temperature given by 

(8) 



'-th ~ /-— 5 



Here i// is the expansion rate during inflation. 

2.2 Scalar field squared expectation value of the inflaton de- 
cay product background 

When discussing the evolution of the r.h. sneutrino condensate, we will need the 
expectation value < >, where ^ represents a generic real MSSM scalar in the 
inflaton decay product background [13]. 

(i) Unthermalised Decay Products during IMD. If we consider the average 
momentum of the scalar modes in the inflaton decay product background to be ~ A;, 
then the average energy density of a real massless scalar held is 

< >=< ^ ^(v*)^ >^k^ <^^> . (9) 

Thus with k !v ~ ms/2 and < p-q, >— f^nPd, where pa ~ {Td/H)ps is the energy 
density of the inflaton decay products during IMD and is the fraction of the total 
inflaton decay product energy density in the real scalar fleld we obtain 

<*^>«^«i™!^. (10) 



(ii) Unthermalised Decay Products during RD. In this case the energy of the 
dominant decay products red-shifts as ~ ^ whilst p^ oc for relativistic 



decay products, with pd{aR) Ps{clr). Thus < >oc Pd/E^ oc oc H. Therefore 
the same relation, Eq. (10), between > and H also holds during RD. 

(iii) Thermalised Decay Products. In this case the energy density of the inflaton 
decay products is p^ — ^ ^gg'*'^ , where for a real scalar field g{T) — 1. Therefore, with 
k lor effectively massless thermalised particles, we find 

<^r2>f«7TT2 ; 7T = ^. (11) 

We note that in the case where the inflaton decay products thermalise immediately at 
the end of inflation, the temperature during IMD is related to H by [20] 

T^K{Mp,HTlYl^ ; (12) 

3 Right-Handed Sneutrino Condensate Evolution 

3.1 Neutrino masses and the r.h. sneutrino scalar potential 

For simplicity we will consider a single neutrino generation. The superpotential of the 
r.h. neutrino superfield, A^, is given by 

W.^X,NH^L+^^ , (13) 

where Hu and L are the MSSM Higgs and charged lepton superfield SU (2)x, doublets 
[2]. The corresponding scalar potential for the r.h. sneutrino is then V{N) = -^-N"^ 
(with N a conventionally normalised real scalar field), where = ml + M% + ml^^. 
Here is the conventional SUSY breaking mass squared term (mo ~ 100 GeV) 
whilst mgjrj. = < > is the effective mass squared term due to the interaction of 
the r.h. sneutrinos with the inflaton decay product backround [13]. In addition, we 
expect terms due to Planck-scale suppressed interactions. For now we will consider 
the evolution of the sneutrino condensate in the absence of such terms. 

If Miv = we will have Dirac neutrino masses, with mv — X^Vu- The r.h. sneutrino 
mass is then given by = + ra'^eff- ^ X^Vu (where v^—K >) we will 
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have Majorana neutrino masses from the see-saw mechanism, my = -j^, such that 
the Yukawa couphng X^, can be expressed as a function of rriy, 

K^(^f. (14) 

The usual idea of the see-saw mechanism [1] is to consider the magnitude of X^, to be 
similar to the charged lepton Yukawa couplings, which requires e.g. Mjv ~ 10^ GeV 
for Xi, ^ Xr ^ 10~^ and m,^ ^ O.leV. However, other natural mass scales could also 
be of interest, for example Mn ~ 100 GeV — 1 TeV, as suggested by the electroweak 
scale and by the scale of the SUSY mass term /iHuHd of the MSSM superpotential [2]. 

3.2 Conditions for r.h. sneutrino to act as a curvaton 

R.h. sneutrino oscillations begin once the r.h. sneutrino mass satisfies ^ H. We 
denote the scale factor at this time by Qosc and the homogeneous sneutrino expectation 
value by A^osc- During inflation, in order to serve as a curvaton, the r.h. sneutrino must 
be effectively massless, <S if/. The quantum fluctuation of an effectively massless 
r.h. sneutrino mode at horizon crossing is given by 5N fa Hi/2ti [20]. In the case 
of a scalar potential consisting purely of a mass term (oc A^^), once the perturbation 
mode is stretched outside the horizon by the expansion of the Universe, its amplitude 
on sub-horizon scales will evolve in the same way as the homogeneous field. This can 
be seen by considering N — No + SN, where Ng is the homogeneous field and SN is a, 
perturbation of wavenumber k. The equations of motion for these are 

No + ^HNo^-V'{No) (15) 

and 

SN + 3H6N - —SN = -V" (N^6N . (16) 

For a mode outside the horizon, the k^/a^ term (<S H'^) will effectively play no role 
in the evolution of the scalar field. Therefore, for V{N) oc A^^, Eq. (15) ^ Eq. (16) 
under ^ 6N. Thus No and SN will evolve in the same way. Therefore SN/N 
for a superhorizon perturbation will be fixed by its value at the onset of oscillations. 



{SN/N)osc- Once coherent oscillations of the r.h. sneutrino begin, the energy density in 
the r.h. sneutrino field will be proportional to its amplitude squared. Therefore if the 
energy density of the Universe becomes dominated by the r.h. sneutrino oscillations 
before the sneutrinos decay, the energy density perturbation when a given mode re- 
enters the horizon will be given by 



(A more precise calculation gives the same result up to a factor of the order of 1 [9].) 
For Hj and Nj constant this corresponds to a scale-invariant perturbation spectrum. 
In order to account for the observed CMB temperature fluctuations, we then require 
that Sp 10-^ [22]. 

3.3 Condensate evolution without Planck-suppressed terms 

The evolution of the r.h. sneutrino expectation value depends on the inflaton decay 
product background and reheating temperature. We will consider the case where 
the r.h. sneutrino begins coherent oscillations during the IMD epoch, aosc < ciR- 
Oscillations begin at Hose ~ '^T'jv ~ 100 GeV. The condition Hr < Hose ~ then 
implies an upper bound on Tr, 



In particular, if Tr ^ 10^ GeV, as would be required by the thermal gravitino upper 
bound [21, 24] if the inflaton decay products are thermalised at the onset of radiation 
domination, then r.h. sneutrino oscillations would generally begin during IMD. How- 
ever, it is possible that thermalisation of the relativistic decay products could occur at 
Tth ^ 10^ GeV even though Tr > 10^ GeV, in which case the thermal gravitino bound 
on Tfi could be evaded. 

We will also assume that thermalisation of the inflaton decay products occurs after 
the Universe becomes radiation dominated. This is true if Eq. (5) is not satisfied, 
which in turn requires that the infiaton mass satisfies 




(17) 




(18) 




(19) 
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where a^c, being due to typical MSSM couplings, is not expected to be much smaller 
than 1. We finally assume that is dominated by the time- independent terms 



N 

when the sneutrino oscillations begin, ^ = ml + M^. This requires that 



^ mj:j . Using < -> unthermahsed background during 

IMD, Eq. (10), this requires that 



m's > ^ . (20) 

For the case of Dirac neutrinos this implies that 

i 1.1 X lOVi'^ (}^^\" f-?}^) (J^) GeV , (21) 

where we have assumed Vu ~ 100 GeV, whilst for the case of Majorana neutrino masses 

This assumption will give the largest r.h. sneutrino energy density at late times for 
a given Tr and N^sc, since the effect of having mg// > mfj^ at H ^ m^^ would be to 
cause r.h. sneutrino oscillations to begin earlier and so to experience a greater dilution 
of the r.h. sneutrino energy density due to expansion. From now on we will consider 

A fundamental condition for the r.h. sneutrino to play the role of a curvaton is 
that the r.h. sneutrinos decay after the Universe becomes dominated by the energy 
density in their coherent oscillations. The r.h. sneutrino decay rate is given by 

The r.h. sneutrinos decay once Fj^^ H. At the onset of r.h. sneutrino oscillations, 
Pjv ~ whilst PS ^ 3H^g^M^ Sm^M^. During IMD, the energy density 

in the infiaton oscillations and the r.h. sneutrino oscillations are both evolve as a~^. 
Thus Pf}/ ps is constant. Once the Universe is radiation dominated, the energy density 
in the dominant relativistic background, p, evolves as whilst the energy density in 
the r.h. sneutrino evolves as a~^. Therefore once a > qr, 

Pn _ f a \ f Pn\ ^ f a \ ^^24) 



anJ \Ps J \aRj 6M'^ 
9 



Thus the Universe becomes r.h. sneutrino dominated once H < Hdom, where 

H^om - (^^J Hn . (25) 

The condition that r.h. sneutrino decay occurs after r.h. sneutrino domination, Fjy^ < 
Hdom, then imphes that 

- [ XlMp, ) [gM^) ■ ^^^^ 

We will refer to Eq. (26) as the late decay condition in the following. 

For the case of Dirac neutrino masses the late decay condition becomes 



which implies that 



0.1eV\V Tp \ViVn.,.^' 



; 5.8 X lO^' [^) i^-^) i^^) GeV . (28) 

Thus in order to have ^ rUo ^ 100 GeV we require that Nosc/M ^ 1 x 10~^. We 
will see that values of Nosc/M in this range require some suppression of Planck-scale 
suppressed non-renormalisable corrections to the r.h. neutrino superpotential. 
For the case of Majorana neutrino masses the late decay condition becomes 

/ m^y /V4^y' Y*2vJiy'' (29) 



such that 



AT/ \ iiLu / \ ^y^pi 



m 



---°=^Se)"^(^)"^(T5^)(^V0eV. (30) 



We will restrict attention to r.h. neutrino masses Mjv ^ mo ^ 100 GeV, since masses 
smaller than this would have little motivation from neutrino mass models or natural 
particle physics scales. In this case m^y fa Mat. Eq. (30) then shows that a right- 
handed Majorana sneutrino curvaton with Mn ^ 100 GeV is possible only if Nose 
is not very small compared with the reduced Planck scale M. {Nosc/M ^ 0.01 if 
Tr'Z 10^ GeV). It also shows that Tr cannot be very small compared with the thermal 
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gravitino upper bound of 10 GeV if Nose ^ M. The requirement that Nose is not very 
small compared with M is a strong constraint on a Majorana r.h. sneutrino curvaton, 
since, as we discuss below, it requires a high degree of suppression of Planck-scale 
suppressed contributions to the r.h. sneutrino superpotential. 

Prom Eq. (17) and the density perturbation constraint, Sp 10~^, the expansion 
rate during inflation is Hi = Tr5pNi > 10^ GeV if Nj fa Nose ^ 1 x 10~^M, correspond- 
ing to the case of Dirac neutrino masses. This can be substantially lower than the 
typical value of the expansion rate during inflation found in conventional SUSY infla- 
tion models, Hj fa 10^^ GeV [25, 26, 27]. Por the case of Majorana neutrino masses 
with Ni Nose ^ O.OIM, the corresponding bound on Hj is Hj > 7 x 10^^ GeV. 
These bounds assume that N does not evolve significantly from the end of inflation to 
the onset of sneutrino oscillations, so that Nose ~ Nj. We will see that the range of 
allowed Hj can be increased if this assumption is altered. 

3.4 Effect of Planck-suppressed terms 

In general, in addition to the globally SUSY scalar potential we expect (in the absence 
of specific symmetries) terms suppressed by powers of the reduced Planck mass, M — 
Mpi/\/Sn, corresponding to the natural scale of supergravity (SUGRA) corrections 
[3, 19]. Thus we expect Planck-scale suppressed non-renormalisable terms to appear 
in the r.h. neutrino superpotential. We also expect contributions to the mass squared 
term of the form cH^, where |c| is model-dependent but expected to be of the order 
of 1 in the simplest models. These arise from terms in the full Lagrangian of the form 

Jd^eS^SN'fN = iFsp/M^, where S is the inflaton field (or any other field with a 
non-zero P-term contributing to the energy density of the Universe) . 
(i) cH^ corrections 

During inflation the value of |c| is constrained by the deviation of the curvaton 
perturbation from scale-invariance. We consider the r.h. sneutrino potential during 
inflation to be 

V{N) ^ -cH^N^ . (31) 
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The expansion rate during inflation is then 

-'-^(--^) 

where ps is the energy density of the inflaton field, which we assume to be constant. 

The index n of the perturbation spectrum as a function of present wavenumber k 
is given by 

2kdSp 

n — 1 

such that Sp/p oc and n — 1 corresponds to scale- invariance. 5p will remain 

constant once the perturbation is outside the horizon, since both N and 5N evolve in 
the same way for V{N) oc N"^. Therefore we have 

where ax is the scale-factor at which a perturbation of wavelength A exits the horizon. 

The value of A?" as a function of the scale factor during inflation is given by the 
solution of 

TV + 3HN = -cH^N (35) 



With H (X a ™, this has a solution of the form N (x a'^, where 

1 



^=2 



(3-m) + J {3 -my -Ac 



(36) 



Thus during inflation (m = 0) the solution corresponds to 

1 



^=2 



-3 + V9^^ ■ (37) 



If |4c| -C 9 then during inflation 7 a; — c/3 and so N/Nx = {a/ ax)^^^^ , where A^^ is 
the value at ax- Thus for a given N and a (for example, their values at the end of 
inflation), Nx oc a^'^^^. The wavenumber at present is related to the scale factor at 
horizon exit (A ~ HJ^) hy k — 2TrHiax/ap, where Op is the scale factor at present. 
Thus 

dNx dNx dax cNx 



dk dax dk 3k 
12 



(38) 



For the case where the energy density is dominated by ps, we have 



Therefore 



^^-JL. (39) 



2k d6p dNx 2c 
n = 1 + — — f — 1 + — . (40 

5, dNx dk 3 ^ ^ 





An 


\c\ < 0.15 


0.1 





Note that c > (< 0) results in a blue (red) spectrum of perturbations. Observation 
requires that |An| < 0.1 [22]. Therefore during inflation we must have 

(41) 

If |c| is close to this upper bound then a significant blue or red perturbation spectrum 
is expected. We note that this would allow the r.h. sneutrino curvaton scenario to be 
consistent with the recent observation by WMAP of a blue perturbation spectrum on 
comovmg scales of the order of 500 Mpc, with n = l.lOlg:^^ [23]. 

Thus |c| is constrained to be not much larger than 0.1 during inflation. This is 
significantly larger than the value |c| ~ 1 expected on dimensional grounds in SUGRA 
models [18]. This problem is similar to the conventional 'ry-problem' encountered in 
SUSY inflation models [25, 26, 27], where order corrections lead both to a large 
deviation from scale-invariance and to insufficient slow- roll inflation [18]. It has a 
natural solution in the case of SUSY D-term hybrid inflation (driven by the energy 
density of a Fayet-Illiopoulos D-term) [25], in which case \F\ =0 and so |c| =0 
during inflation (although a non-zero value is expected once inflation ends and coherent 
inflaton oscillations begin, since Jd'^OS^S^^^ = (l^;.'^!^ + |-^5p + -) |*|^ with $ a 
general scalar superfield). Since in the r.h. sneutrino curvaton scenario inflation must 
still be driven by an inflaton, it is possible that the mechanism which suppresses |c| 
for the inflaton during inflation also suppresses |c| for the curvaton. 

After inflation we generally expect a non-zero F-term from the energy density of the 
coherently oscillating inflaton fleld. The solution of Eq. (35) during IMD corresponds 
to Eq. (36) with m — 3/2. Thus assuming that N has a constant value during inflation, 
Ni, we flnd 

Nose faoscV f Hi \'"'/' ^^2) 



"osc 

Nt \ ttp J \Hosc 
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where Hose ^ rrij^. 

We previously derived lower bounds on the value of Hi compatible with the ob- 
served density perturbations based on the assumption that Nose ~ Nj. However, if 
c < during IMD then the growth of N after inflation will allow a wider range of Hj 
to be compatible with a given value of Nose, 

' nSpNosc^ 



27 



In 



Hi 



3 In(^) ' 



(43) 



where we have used Nj ^ Hj/nSp in Eq. (42). For example, for the case of Dirac 
neutrino masses, with 5p fti 10~^, 100 GeV and Nose ~ 10~^M, it is possible to 

have Hj 10^ GeV if 7 = 1.45, corresponding to c = —4.3. This also shows that |c| 
need not be small compared with 1 after inflation. 

We conclude that the cH'^ correction to the r.h. sneutrino mass squared term is 
no more problematical for the curvaton than for the conventional inflaton: we require 
that |c| ^ 0.1 during inflation, in which case a signiflcantly blue or red perturbation 
spectrum can arise depending on the sign and magnitude of c. After inflation |c| 
need not be small compared with 1 (in the case where c < 0) and the cH^ term can 
even widen the range of expansion rate and energy density during inflation which is 
compatible with the observed density perturbations, 
(ii) Non-renormalisable superpotential corrections 

We next consider adding Planck scale-suppressed non-renormalisable terms (NRTs) 
to the r.h. sneutrino superpotential, 

^" = ^^"^^ + ;si^- 

where n\ is a symmetry factor and Xn 1. For large enough \N\ the scalar potential 
becomes dominated by the NRT contribution. 

The effect of the NRTs is to place an upper limit on the value of Nose- During inflation, 
the effective mass V" (N) must be small compared with H^. After inflation and during 
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IMD, H'^ oc a^^. Therefore ^"(A'^) may become larger than if^, at which point the 
N field will begin to evolve. will then track the value at which V" {N) ^ iJ^, since 
as decreases below V" {N)^ the rate of roll of N will increase until the rate of 
decrease of V" {N) matches rate of the decrease of H'^. In particular, when N will first 
begins to slow- roll we have ZHN —V'{N). For V{N) oc N'^^'^~^\ this has a solution 
N oc a^/^'^~'^'^\ such that V" {N) oc oc if^. So the slow-roUing solution will be such 
that V" {N) tracks and is of the same order of magnitude as H^. Thus the value 
of N at which V" {N) a; H"^ places an upper limit on N for a given value of iJ, N^^^ 
given by 



where is a constant of order 1. Therefore for Nu^ to be greater than Nose at 
we require that n is greater than n^j^, where 

In i ^) 



For the case of Dirac neutrino masses, we require from the late decay condition 
that Nosc/M > 1 X 10"^. Thus with Hose ^ rrij^ ^ 100 GeV, Nosc/M ^ 1 x 10'^ 
and An ~ q;„ ~ 1, we require that n > num — 5.3. Thus a suppression of Planck- 
scale suppressed NRTs in the r.h. neutrino superpotential of dimension less than 6 is 
required in this case. This might be achieved by a modest discrete symmetry. 

For the case of Majorana neutrino masses, with Hose ^ = 100 GeV, Nosc/M pa 
0.01 and A„ w a„ ~ 1, we require that n > = 10.2. This is a significant problem 
for a Majorana r.h. sneutrino curvaton. It requires a high degree of suppression of 
Planck-suppressed non-renormalisable terms, eliminating all NRTs in the r.h. neutrino 
superpotential up to n < 11. However, as we will discuss, there is likely to be a more 
severe problem for the Majorana r.h. sneutrino curvaton, namely the survival of the 
r.h. sneutrino condensate in the inflaton decay product background. 

Although Planck scale-suppressed NRTs appear to disfavour a r.h. sneutrino cur- 
vaton, it should be noted that conventional SUSY inflation models also have problems 




(46) 




(47) 
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with Planck scale-suppressed NRTs. Chaotic inflation models require that the infla- 
ton expectation value is greater than M [27], whilst SUSY hybrid inflation models 
require that N is close to M when scales corresponding to observed cosmic microwave 
background (CMB) perturbations exit the horizon (assuming natural values of the 
renormahsable gauge and Yukawa couplings) [25, 26]. In this sense the r.h. sneutrino 
curvaton may be no more problematical that conventional SUSY inflation models with 
respect to Planck scale-suppressed superpotential terms. 

3.5 Decay temperature of the Dirac and Majorana r.h. sneu- 
trino condensate 

The neutrino mass is related to the Yukawa coupling by rrii, — XuVu- The r.h. sneutrino 
mass is simply given by the SUSY breaking mass term, mf, — ruo ^ 100 GeV. Thus 
the temperature of the Universe when the condensate decays, Tjy^, is given by ^ 
H{Tfj^), where H{T) = kxT^Mpi and kr = {47r^g{T)/45y/\ Thus with the r.h. 
sneutrino decay rate given by Eq. (23), the temperature of r.h. sneutrino decay is, in 
general, 

- [ inkr-, J ■ ^^^^ 
For the case of the Dirac r.h. sneutrinos, Eq. (48) imphes that 



where we have used v„ ~ 100 GeV and /c^^^ ~ 5, corresponding to 7, and Vi 
{i = 1,2,3) as light degrees of freedom [20]. Thus the Dirac r.h. sneutrino condensate 
typically decays in the temperature range 1 — 10 MeV. (Note that if the energy density 
of the r.h. sneutrino condensate dominates the energy density of the Universe when 
it decays, T^^ should then be interpreted as the temperature to which the Universe 
reheats after the condensate decays.) Since the Dirac r.h. sneutrino condensate de- 
cays well below the temperature at which weakly interacting particles of mass of the 
order of rriw freeze out of chemical equilibrium {Tfreeze ~ 1 — 10 GeV), an important 
constraint on the Dirac r.h. sneutrino curvaton is the requirement that the lightest 
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supersymmetric particles (LSPs) produced in decay of the r.h. sneutrino condensate 
do not overdose the Universe. This requires that the LSPs decay before nucleosyn- 
thesis at T ^ 1 MeV, in order that the light element abundances arc not disrupted by 
photo-dissociation due to LSP decay cascades [24]. Therefore the LSP lifetime must 
satisfy 

tlsp t H-\T ^ 1 MeV) ^ Is . (50) 

Thus although the LSP must be unstable its lifetime can be much longer than the time 
required to escape particle detectors (~ 10~^s), in which case experimental searches 
for SUSY particles would be unaffected. A wider range of LSP candidates would be 
allowed than in the R-parity conserving case where the LSP properties are constrained 
by the thermal relic cold dark matter density. This might be testable at future colliders 
such as the CERN Large Hadron CoUider. 

For the case of Majorana r.h. sneutrinos, Eq. (48) implies that 

^..-f^^V' = 2.2xI0'r^V'7^) GeV. (51) 



47rvf,kT- y ■ VO.l eVy VlOO GeV 

Thus since m ^ typically will not be much larger than 1 TeV as a result of the late decay 
condition, Eq. (30), we expect that T^^ ~ 10"^ GeV for the Majorana r.h. sneutrino 
condensate. 

3.6 Thermalisation of the right-handed sneutrino condensate 

We have been assuming throughout the preceeding discussion that the r.h. sneutrino 
condensate survives in the enviroment of the inflaton decay products. However, it 
is possible that r.h. sneutrinos in the condensate could be inelastically scattered by 
collisions with the MSSM particles in the inflaton decay product background. If the 
relativistic particles in the background have energy E^, then the scattering cross- 
section of a r.h. sneutrino, at rest in the condensate, from a relativistic particle in the 
inflaton decay product background via Higgs(ino) or (s)lepton exchange is expected 
to be 
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where a,, — A^/47r, ag — j^'K-, where is a typical MSSM gauge/Yukawa couphng 



and EcM ~ \J Edrrifj is the centre of mass energy of the process. The scattering rate 
from the decay product background is then F^^^ = na^j^^, where n is the number 
density of particles in the inflaton decay product background. The largest scattering 
rate will occur for the largest n and smallest E^, which corresponds to thermahsed 
decay products (T < Tth). In this case n = g{T)T^/7r^ and E^^T [20]. Therefore 

_ g{T)a,ag 

^Nsc ^^3) 

Then assuming that thermalisation of the inflaton decay products occurs during RD, 
the condition that the condensate is unthermalised is F < H{T) = hrT^Mpi, which 
implies 

A.<^^V'V^V^ (54) 



agg{T)J \Mpi 
For the case of Dirac neutrino masses, this condition becomes 

~ hrr^ ~ ^ -.nn TT^^TTTTTT ^eV, (55) 



agg{T) J \MpiJ q;^^ VlOO GeV/ VlOO GeV. 

where we have used g{T) ^ 200 and /cj^ 20. Therefore, since m,^ ~ O.leV, the Dirac 
r.h. sneutrino condensate will not be thermahsed. 

For the case of Majorana neutrino masses, the non-thermalisation condition be- 
comes 



< Jl . 121! f^^VaV . (56) 



agg{T) J Mpi ' ag VlOO GeV 
Thus typically we expect that this will not be satisfied. For example, if we consider 
ag to correspond to the top quark Yukawa (r.h. sneutrinos scattering from thermal 
top quarks via H,^ exchange) then ^ 0.1 and so uiy ^ 10^'^cV would be required 
to evade thermalisation. Therefore with rriy ~ O.leV the Majorana r.h. sneutrino 
condensate will be thermahsed at T^h, when the inflaton decay product background 
thermalises. 

One possible escape from this conclusion is that the inflaton decay product back- 
ground could remain unthermalised until the energy density of the sneutrino con- 
densate comes to dominate the energy density of the Universe. Since the Majo- 
rana r.h. sneutrino condensate decays at a temperature typically around 10^ GeV, 
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if Tth ^ 10^ GeV then the Majorana r.h. sneutrino curvaton may remain a possibihty. 
For the case of perturbative inflaton decays, the thermahsation temperature of the in- 
flaton decay product background Tth, Eq. (8), is proportional to (Tji/ms)^. Therefore 
a very low background thermahsation temperature is a possibility if Tr <^ ms- 

The above applies to condensate thermahsation during RD. It is straightforward to 
see that if condensate thermahsation does not occur during RD then it will not occur 
earlier during IMD. If we assume that the background is thermalised during IMD 
(which should give the largest rate of condensate thermahsation), then since during 
IMD r^gp oc and H (xT^,it follows that once T > Tr the thermahsation condition 
^Nsc ~ ^ '^i^^ become more difficult to satisfy. Thus if the r.h. sneutrino condensate 
does not thermalise during RD it will not thermalise at all. 

In discussing the thermalisation rate, we have implicitly assumed that the mass 
of the and L fields coupling directly to the r.h. sneutrino condensate is small 
compared with T and Ecu ~ {Tmffy^'^. The efi^ective mass is given by X,^ < N >, 
where < > is the amphtude of the coherent oscillations. In general < > is given 

Thus for the case of Majorana r.h. sneutrinos the effective mass is given by 
which implies that 

x-^(a^)"(To4vr 

where we have used kx ~ 20 and f „ = 100 GeV. Thus for Nggc/M < 1, ^ 100 GeV 
and for a decay temperature of 10^ GeV for the Majorana r.h. sneutrino, we find that 
A,^ < > is smaller than T and Ecm and so may be neglected when considering 
condensate thermalisation. 
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For the case of Dirac r.h. sneutrinos the effective mass is given by 



which imphes that 



X^<N >=3x 10 
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0.1 evy Vio Mevy 



f Tr /lOOGeVX fNosA ^ 

Thus the effective mass is generally negligible in the case of Dirac r.h. sneutrinos. 



4 Conclusions 

We have considered the possibihty that a r.h. sneutrino could play the role of a 
curvaton in the cosmology of the MSSM extended to accomodate neutrino masses. 

In the case of a Dirac r.h. sneutrino, the expectation value of the r.h. sneutrino 
at the onset of its coherent oscillations must satisfy Ngsc/M ^ 10"^, in order that 
the energy density of the r.h. sneutrino condensate dominates the Universe when it 
decays. As a result, Planck-scale corrections to the r.h. neutrino superpotential must 
be suppressed, eliminating all superpotential terms oc A^" with n < 6. cH^ corrections 
to the r.h. sneutrino mass squared must also be suppressed during inflation (|c| ^ 0.1). 
The inflaton sector of the Dirac r.h. sneutrino curvaton scenario can have a much 
smaller expansion rate during inflation than conventional SUSY inflation models, with 
Hi 10^ GeV possible. (Including the effect of a negative cH"^ term after inflation, 
values of Hj as small as 10® GeV or less arc possible.) 

In addition, depending on the sign and magnitude of the cH^ correction, it is 
possible to have a signiflcantly blue (c > 0) or red (c < 0) perturbation spectrum. The 
recent suggestion from WMAP observations [23] of a blue perturbation spectrum on 
comoving scales of the order of 500 Mpc could therefore be accomodated within the 
curvaton scenario. 

The late decay of the Dirac r.h. sneutrino condensate (at T fti 1 — 10 MeV) requires 
that R-parity be broken and that the LSP decays before nucleosynthesis, corresponding 
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to a lifetime shorter than Is. (Since this would result in the loss of LSP cold dark 
matter, a new dark matter candidate would also be needed.) Thus the Dirac curvaton 
scenario predicts that LSP properties will typically be inconsistent with thermal relic 
cold dark matter. 

In the case of a Major ana r.h. sneutrino curvaton, we find that the requirement that 
the r.h. sneutrino condensate dominates the energy density of the Universe before it 
decays implies that Nosc/M > 0.01. This imposes a strong constraint on Planck-scale 
suppressed contributions to the r.h. neutrino superpotential, requiring elimination 
of all terms oc with n < 11. However, a more severe problem may arise from 
scattering of the condensate sneutrinos by particles in the thermal background, since 
it is likely that the Majorana r.h. sneutrino condensate will be thermalised as soon 
as the inflaton decay product background thermalises. Thus if the inflaton decay 
products thermalise earlier than the time of r.h. sneutrino condensate domination of 
the energy density then the Majorana curvaton scenario will be ruled out. Only a 
sufficiently low background thermalisation temperature could evade this conclusion. 
Therefore, although not absolutely excluded, the Majorana r.h. sneutrino curvaton 
scenario appears to be more difficult to implement that the Dirac curvaton scenario. 

In conclusion, wc find that a r.h. sneutrino can serve as the source of the density 
perturbations leading to structure formation. For the favoured case of a Dirac r.h. 
sneutrino, this requires some suppression of its Planck-scale suppressed superpotential 
self-interactions together with sufficiently rapid R-parity violation. If these conditions 
are met, the Dirac r.h. sneutrino would provide us with the only curvaton candidate 
which is strongly motivated by particle physics. 

Note Added: After completing this work we became aware of [28] and [29], which 
also discuss the case of a Majorana r.h. sneutrino curvaton. 
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